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We present a truncation scheme of the effective average action approach of the nonperturbative 
renormalization group which allows for an accurate description of the critical regime as well as of 
correlation functions at finite momenta. The truncation is a natural modification of the standard 
derivative expansion which includes both all local correlations and two-point and four-point irre- 
ducible correlations to all orders in the derivatives. We discuss schemes for both the symmetric and 
the symmetry broken phase of the O(N) model and present results for D = 3. All approximations 
are done directly in the effective average action rather than in the flow equations of irreducible 
vertices. The approach is numerically relatively easy to implement and yields good results for all N 
both for the critical exponents as well as for the momentum dependence of the two-point function. 

PACS numbers: 05.10.Cc, ll.10.Gh, 64.60.ae, 64.60.fd 



The nonperturbative renormalization group (NPRG) 
technique is based on an exact flow equation of the ef- 
fective average action (or generating functional of irre- 
ducible vertices) 0, 0] and has been applied to a large 
variety of systems, see e.g. Refs. [3|-|6( for reviews. It 
proved especially useful when applied to critical phe- 
nomena where often even relatively simple truncation 
schemes yield an accurate description of the critical re- 
gion [H, 0| . While the exact flow equation of the effective 
average action can almost never be solved, it allows for 
novel nonperturbative approximation techniques. One 
successful approximation strategy is the derivative ex- 
pansion, where the effective average action is expanded 
consistently to a given order in spatial derivatives, but 
no truncation is made in the power of the fields. The 
derivative expansion has been applied with success to 
O (TV) models 

The derivative expansion allows strict control over the 
symmetries of the studied models since all approxima- 
tions are done in the effective action which is expanded 
in invariants of the model. This approach automatically 
yields flow equations which are both closed and further 
obey the symmetry of the original model. In contrast, if 
a direct field expansion of the effective action is employed 
and approximations are done at the level of the flow equa- 
tions of irreducible vertices, the invariance of the action 
is not guaranteed. In general, approximations at the level 
of flow equations of vertices are therefore more difficult 
to control. On the other hand, the derivative expansion 
can only access the asymptotic small momentum regime 
of the theory and it cannot be applied to calculate cor- 
relation functions at finite momenta. 

Recently, several approximation schemes which allow 
to calculate correlation functions at finite momenta were 
developed [HT^. The most sophistacted of these is the 
scheme presented in Refs. [l(| El (BMW scheme) which 
is based on an approximate solution of the exact flow 
equation of the two-point vertex in presence of a back- 
ground field. All of these approaches rely however on ap- 
proximations at the level of flow equations for irreducible 
vertices. Instead, here we want to develop a scheme in 



which all approximations are done directly at the level of 
the effective average action. This allows for a transpar- 
ent calculation of the momentum dependence where the 
full symmetry of the model is always obeyed by the flow 
equations. While controlled, such a scheme does how- 
ever not neccessarily guarantee more accurate results for 
the critical exponents. The scheme we discuss below is 
a natural modification of the usual derivative expansion 
and is based on a local potential which is supplemented 
by a momentum dependent potential which accounts for 
nonlocal correlations up to the four-point vertex. Sim- 
ilar schemes were previously used to calculate the one- 
particle spectral function of Bose condensates [H, 03] 
(see also [H EH) and also the thermal fluctuations of 
crystalline membranes such as graphene [l7j (see also 
[la])- Here, we extend the scheme to include, besides 
the nonlocal terms, the full local potential and test it on 
the O(N) model, see e. g. Ref. [19| for a recent summary 
of results on the O(N) model. We develop two schemes, 
a nonlocal potential approximation (NLPA) for the or- 
dered state and a NLPA for the symmetric state. Both 
schemes allow for an investigation of the critical region. 
We first introduce in Sec. fl] the NLPA approach for both 
the ordered state and the symmetric state. We present 
results from a numerical solution of the flow equations in 
Sec. HIl with results for the critical exponents presented 
in subsection III Al In subsection III Bl we assess the qual- 
ity of the approach in the finite momentum regime. We 
conclude in Sec. IIIII 



NONLOCAL POTENTIAL APPROXIMATION 



We begin with an approximation scheme for the or- 
dered state, where we explicitly incorporate a finite order 
parameter into the invariant effective action. 
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A. NLPA for the ordered state 

Starting point of the NLPA approach is an effective av- 
erage action which consists of both a nonlocal potential 
term, which is restricted to second order in the invariant 
densities and is characterized by the coupling function 
u A (k), and a local potential term U A (p — p A ) which may 
be an arbitrary function of p — p A , where p A is the (cut- 
off dependent) order parameter density and p = tp 2 /2, 
where cp is a field with A^-components. To avoid double 
counting of correlations, we define u A (k) to be completely 
nonlocal with ma(0) = 0. We furthermore keep the full 
momentum dependence a A (k) of the quadratic term in 
the action, and thus approximate the effective average 
action, after subtraction of the nonintcracting contribu- 
tion (1/2) J k G^ A (k)(p k ■ (f-k, as 

TaM = \ j [ cr A( fc )Vfc ■ V-k + u A (k)Ap k Ap- k ] 

+ [u A (p-p° A ), (1) 



where Ap*. is the Fourier transform of p(x) — p A . We 
use the notation J k = J d D k/(2ir) D and J — J dPx 
for integrals over momenta and integrals over coordinate 
space, respectively. U a (t) can, for finite cutoff A, be 
expanded in r, 



(n). 



(2) 



with r = p — p A or, for the symmetric scheme discussed 
in subsection II Bl r = p. Note that the effective action 
(TTJ) docs not contain all terms of a complete derivative 
approximation to order q 2 , in which the derivative term 
of the action would also include an expansion to all pow- 
ers of p — p A . In the present scheme, one could easily 
improve upon the action (TTJ) by including e.g. also ad- 
ditional terms which can be parametrized by only one 
momentum. One obvious extension would be to include 
a term J xy [p{x) - p° A } 2 [p(y) - p° A }n{x - y), with some 
function k(x) which would be a generalization of the 
(d l _ l p) 2 (p — p A ) term encountered in a derivative expan- 
sion. Such an extension is both straightforward and nu- 
merically feasible. 

At the same time, the present approach goes well be- 
yond the derivative expansion in that it includes the full 
momentum dependence in the first two terms of Eq. (TTJ) . 
As in the derivative expansion, the effective average ac- 
tion obeys the full O(N) invariance throughout the en- 
tire flow. To determine the flow of U A , we can use the 
standard technique [11] and evaluate the flow of T A [<p] 
for a homogeneous (a;- independent) field <p such that 
V- l T A [(p] = U A (p - p%) with p = ip 2 /2 and where V 
is the volume. We shall now assume that N > 2, so that 
there is at least one gapless transverse mode. The flow 



of the local potential is then given by [l[ 

d A U A {p- P l) = \ ( d A R A (k){G Atll (k,p) 



where 



+ (N-l)G A>x (k,p)}, (3) 



G- K \{k,p) = a A (k) + U' A (p - p° A ) + G X(k) , (4a) 
GXJ(M = <7 A (*) + 2p[u A (k) + UUp-pU 

+ U' A (p- P ° A ) + G a X(k). (4b) 

Here, the cutoff regulated noninteracting Green's func- 
tion is 



R\(k) . 



(5) 



and R A (k) is a regulator for small momenta with k < A. 
The only difference of Eq. Q to the standard form used 
in a derivative expansion of T A is the presence of the full 
functions u A (k) and a A (k) in Eqs. (|4"a|) and (|4bl) rather 
than just their leading terms of a fc-expansion. To de- 
termine the flows of cr A (k) and u A (k) we invoke a field 
expansion of Ta[^] in terms of A.ip^ = (pf, — tp A 6 a iSk,o 
with p A = (ip A ) 2 /2. Here we have, without loss of gen- 
erality, assumed an order parameter field <p A which is 
directed in the a — 1 direction of the internal space. To 
determine the flows of a A and u Al we need the lowest 
order irreducible vertices (up to four-point), which have 
the form 



r A 2 L( fc ; = S ab a A {k) + 2S al S bl p° A u A (k) , (6a) 

r A 3 Lc( fe l' fe 2,fc3) = VA 5 a i5 bc U A {ki) + 5 bl SacU A {k 2 ) 

+ 5 cl S ab u A (k 3 ) + ((p° A ) 3 U A 3) 5 al 6 bl S cl 

(6b) 



r A 4) abc<i( fc l ' • • fe <0 = S ab S c dU A (ki 2 ) + 5acSbdU A (k 13 ) 

+ 5ad5bcU A (k 1A ) + 2p° A U A 3) 5 ab 6 cd 
x (Sai + $ci) + 5 ac 8 bd (5 al + S bl ) 

+ S a dS bc (Sal + Sbl) 

+ 4(p° A ) 2 U ( A ) 6 al 5 bl 5 cl 5 dl , (6c) 

where we defined u A (k) = u A (k) + U A and = \ki + 
kj | . The flow of the order parameter follows from the 



requirement that d A T A 

dAP° A = ' 



(1) 



0. This yields [2C 



2u A (0) 



U A (0) + 2u A (q) + 2p° A U ( A 3) ]G A ^(q) 



+ (N-l)u A (0)G A , x (q)}, (7) 



where G A ^ a (k) = G A . a (k, p A ) for a =_L, || and G A ^ a (k) = 
—G Aa (k)d A R A (k). The flow of a A (k) follows from the 
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(2) (21 

flow of T A j_(k) = T A aa (k, —k) where a ^ 1 is a direction 
transverse to the order parameter field. 

d A a A {k) = { {G±(q)u A (q') - G||(g)u A (?)} 

J q 

-2p° A [ \G\\{ q ')G^q)ul{q') 



+ G x (q')G [l (q)ul(q)}, (8) 



and we defined q' = \k + q\. The flow equation of u A (k) 

(2) (2) 

can be obtained from the flow of T A ^ (k) = T A n (fc, — k) 
which reads 

9A^l\(k) = \ I {(N-l)G A , x (q)[u A (0) + 2p° A U^] 

J q 

+ G Ai ||(g)[u A (0)+2«A(9 / ) 

+ 12p A [/f +4(p A ) 2 f/f]} 

-2p° / f(JV-l)6 Ai± (g / )G A) ±(ff)«A(*) 



+ G Aill {q')G Aill (q)[u A {q) + u A (q') 

+ u A (k) + 2 PA U A 3) } 2 } 

+ [« A (0) + 2w A (fc) + 2p A C/| 3) ]a A ^ • (9) 

Combining Eqs. (0), (jHJ, ©, and keeping in mind that 
T^Jk) = a A {k) + 2p° A u A (k), one finds 

d A u A (k) = J [G At \\(q) - G A .±(q)]u A (q') 
G A ,ll(g){2[/f+p A [/| 4) 
-U A 3) [u A (q) + p° A U A ^]/u A (0)} 
' (N - l)G A , ± (q')G A , ± (q)u 2 A (k) 

+ G AA{ (q')G At{l (q) [u A (q) + u A (q r ) 
+ u A (k)+2p A U A 3) ] 2 } 

+ ^{G l{ (q')G x (q)ul(q>) 

+ G x (q')G l] (q)ul(q)]. (10) 

This completes the derivation of the flow equations, 
which are uniquely determined by the effective action 
(JTJ). The flow Eqs. PI7I8I10I) form a closed set which can 
be used to calculate the full momentum dependence of 
the self-energies in a controlled manner and the only ap- 
proximation is the form of the effective action as stated 
in Eq. ([I]). By construction, the approach reproduces 
exactly the correct structure of the leading order per- 
turbation theory, which is dominant at large momenta. 
Also by construction, it reproduces the leading terms in 



a derivative expansion of both u A (k) and a A (k) to lowest 
order in the fields, which dominate the behavior in the 
infrared. The same is true also for the symmetric scheme 
which we discuss below. 




FIG. 1: Dependence of the anomalous dimension r\ on the 
order n of the polynomial approximation of the local poten- 
tial Ua{v) = 2~2 n =o tf /j! f° r the symmetry broken phase 
(upper curve) and the symmetric scheme (lower curve) in 
which the fixed point is approached from within the sym- 
metric phase. Values shown are for N = 2 and D — 3 . 



B. NLPA for the symmetric state 

We now derive flow equations which are valid for the 
symmetric phase, which are even simpler. In the NLPA 
for the symmetric state the distance to the critical point 
is controlled by a mass term r A in the propagator which 
vanishes at criticality in the limit A — > 0. We write the 
Ansatz for T A in the NLPA as 

TaM = ~ J {[ (T A(fc) + r A ]ip k ■ <P-k + u A (k)pkp-k} 



U A {p) , 



(11) 



where pt is the Fourier transform of p(x) = cp 2 (x)/2 

and where we put U A = to avoid double counting 
of the mass term which is already accounted for by r A . 
The action (jlip yields again unique flow equations for r A 
and the functions U A (p), a A {k) and u A {k) which can be 
easily derived. We define the vertices now as expansion 
coefficients of T A around cp — 0. The flow for the two- 

(2) 

point vertex T A ' ab (k, —k) — 6 a f ) 'E A (k) is then 

d A Z A (k) = i J G A (q) [2u A (q') + Nu A (0)] , (12) 

where S A (fc) = r A + u A {k) with a A (0) = and where 
G A (q) = -G 2 A (q)d A R A (q) withC^ = G^ A +£ A (fc). The 
flow of the two-point vertex is in the symmetric phase not 
sufficient to extract also the flow of u A (k) = u{k) + 
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TABLE I: Values for the anomalous dimension 77 for various N and D = 3 from different approaches. The columns correspond 
to the symmetric NLPA, the NLPA for the ordered phase, results from the background field scheme (BMW) 11], the first order 
derivative expansion (DE), field theory (FT), variational perturbation theory (VPT) [2l|] and Monte Carlo (MC). 



N 


sym. NLPA 


ord. NLPA 


BMW 


DE 


FT 


VPT 


MC 





0.042 




0.034 


0.039 [22] 


0.0272(3) [23J 


0.031(1) 


0.0303(3) [24] 


1 


0.042 




0.039 


0.0443 [7] 


0.0318(3) [23] 


0.034(7) 


0.03627(1) [25] 


2 


0.041(5) 


0.049 


0.041 


0.049 [22] 


0.0334(2) [23] 


0.035(6) 


0.0381(2) [26] 


3 


0.040 


0.046 


0.040 


0.049 [22] 


0.0333(3) [23J 


0.035(0) 


0.0375(5) [27] 


4 


0.038 


0.042 


0.038 


0.047 [22] 


0.0350(45) [28] 


0.031 


0.0365(10) [29] 


10 


0.026 


0.024(5) 


0.022 


0.028 [22] 


0.024 [31] 


0.0216 





and we must extract its flow from the four-point vertex. 
This yields 

dAUA(k) = J^G A (q)U A 3) - J^G A (q)G A (q'){ 

(N - l)u A (k) 2 + [u A (k) + u A (q') + u A {q)] 2 } 

G A (q)G A (q){ 

[u A (q') - u A ( g )][fi A (0) + 2u A (q)]} , (13) 

with q' = \k + q\. The flow of the local potential U A (p) 
is given by Eq. ([3]) with p\ = and with a A (k) replaced 
by cr A (k) + r A in Eqs. (|Ia} and (jib]) . 

II. RESULTS 

We have solved the flow equations both in the sym- 
metric phase and the symmetry broken phase numeri- 
cally for D = 3 and different values of N. For D = 3 the 
field expansion of the local potential actually converges 
relatively fast [7| so that one can work with a finite order 
approximation of the local potential. We have used an 
expansion of U A (p) up to eighth order in p in both the 
symmetric and the symmetry broken scheme and have 
checked that the values of the anomalous dimension 77 
are already converged at this level of truncation. The 
convergence can clearly be seen in Fig. [T] where we show 
the values of 77 for different maximal powers of p. All re- 
sults presented below were calculated with all terms up 
to order p 8 . 

For numerical stability we choose an exponential cut- 
off, 



RA(q 2 ) 



exp((7 2 /A 2 



1 



(14) 



where Z^ 1 = l + dk2V A (k)\k=o is the wavefunction renor- 
malization. Usually the prefactor a is tuned in such 
a way as to extremize the critical exponents, e.g. the 
anomalous dimension 



7/ = Ad A In Z A . 



(15) 



This ensures a minimal sensitivity of the results to small 
variations in a Q. In the present scheme we do not 
observe an extremal value of 77 as a function of a. Instead, 
we observe a steady decrease of 77 when a is increased and 
a minimum which is only reached asymptotically for large 
a. For the symmetric scheme, the dependence of 77 on a 
is already essentially flat for a > 5 and we choose a = 5 
for our analysis below. Similarly, in the symmetry broken 
phase only a small decrease of 77 is detected on increasing 
a from 1 to 2 and 77 is then essentially unchanged up to 
a = 3. We fixed a = 2 for the analysis below. 



A. Critical exponents 77 and v 

The value of 77 can be easily determined from the flow of 
the quantity o A (k) and its low momentum structure via 
Eq. (|15l) . To determine the thermal exponent v, we use in 
the symmetric phase the value of the fully renormalized 
mass term r* = limA->o r A which scales as r* ~ (r Ao — 
re) 21 ' where r c is the critical value of the mass term at the 
initial cutoff scale Ao. Similarly, in the symmetry broken 
phase we analyse the scaling of the order parameter = 
lim A _).o p A which scales as p* ~ (pa — Pc) 2/3 where /3 
is the critical exponent of the order parameter and p c 
the critical value of p A at the initial scale Ao. From 
/3 and 77 we can extract v via the hyperscaling relation 
v = 2/3/(L>-2 + 77). 

In Tables U and [IT] we show our results for the criti- 
cal exponents 77 and v and compare them with results 
from various other approaches. Somewhat surprisingly, 
and in contrast to what is observed in a standard deriva- 
tive expansion, the results for the critical exponents are 
generally better in the scheme where one approaches the 
critical point from the symmetric side, where a field ex- 
pansion around p = rather than around a finite value 
p° is employed. As can be seen from Table HI the values 
for 77 in the symmetric scheme are, except for N = 0, 
quite close to those obtained within the BMW scheme 
of Ref. [HI- For large TV, it is known that 77 behaves 
as 77 = 0.21 /N [30], and the result from the symmetric 
scheme for N = 10 is already close to this value. 

The results for the approach from the symmetry bro- 
ken phase are similar to those of the leading order deriva- 
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TABLE II: Values for the anomalous dimension v for various N and D = 3 from different approaches. The columns correspond 
to the symmetric NLPA, the NLPA for the ordered state, results from the background field scheme (BMW) [TlJ, the first order 
derivative expansion (DE), field theory (FR), variational perturbation theory (VPT) [2lJ and Monte Carlo (MC). 



N 


sym. NLPA 


ord. NLPA 


BMW 


DE 


FT 


VPT 


MC 





0.58 




0.589 


0.590 [22] 


0.5886(3) [23] 


0.5883 


0.5872(5) [32] 


1 


0.62 




0.632 


0.6307 [7] 


0.6306(5) [23] 


0.6305 


0.63002(10) [25] 


2 


0.66 


0.68 


0.674 


0.666 [22] 


0.6700(6) [23J 


0.6710 


0.6717(1) [26] 


3 


0.70 


0.72 


0.715 


0.704 [22] 


0.7060(7) [23J 


0.7075 


0.7112(5) [27] 


4 


0.74 


0.76 


0.754 


0.739 [22] 


0.741(6) [28] 


0.737 


0.749(2) [29J 


10 


0.89 


0.89 


0.889 


0.859 [22] 


0.859 [31] 


0.866 





tive expansion (where terms up to order 0{q 2 ) are kept), 
see Table HI A possible reason for the inferior accuracy 
of the scheme for the symmetry broken phase compared 
with the accuracy of the symmetric scheme is that all 
nonlocal correlations are determined already from the 
two-point function whereas in the symmetric scheme the 
nonlocal potential flow is determined from the four-point 
function. Including further terms in the effective action 
is expected to improve also the results of the symmetry 
broken NLPA. 

The results for the thermal exponent v are similar in 
both schemes and generally close to the most accurate 
MC results with deviations never more than about 3%. 
Our values are also close to values from other approaches. 

B. Beyond the universal regime 

Both the schemes for the symmetric and the symme- 
try broken phase reproduce the logarithmic behavior of 
the self-energy at large momenta £(fc) ~ u\ Q \n(k]u\ ) 
which can be derived from perturbation theory [33j. To 
assess the accuracy of the calculated self-energy over the 
whole momentum regime a useful quantity is the small 
u\ limit of the one-dimensional integral (£(z) is the Rie- 
mann zeta function) 

onuAg J q z + T,(q) 

where S(q) is the full self-energy at criticality, 

(2) 

r A,ab( fc : ~ k ) = ^a&£ A (fc) and S(fc) = lim A ^o S A (fc). The 
quantity c is finite in the limit ma /Ao — > and has 
physical significance for N — 2 where it relates to the 
suppression of the critical temperature of the weakly in- 
teracting Bose gas in D = 3 dimensions [33]. The inte- 
gral in Eq. (|16[) is dominated by contributions from the 
crossover regime k ~ u\ where the momentum depen- 
dence of the self-energy changes from the perturbative 
ln(fc) behavior at large momenta to the anomalous fc 2- '' 
scaling at small momenta. The value of c has been es- 
timated for different N from Monte Carlo simulations 
[35M37I ] and has been used to quantify the accuracy of 
various approaches [H, [HI OS} ■ To calculate c we used a 
small initial value of u\ , ua /Aq = 0.001. Again we find 



generally better values for the symmetric scheme. For 
N = 2, the value is about 15% too high when compared 
to MC results and for N = 1 the difference is slightly 
larger. For N = 4 the difference is less than 5%. In com- 
parison with the BMW scheme the differences are 8% for 
N = 2 and rapidly decrease for larger N, see Table Hill 



III. CONCLUSIONS 

We have presented a straightforward nonlocal poten- 
tial approximation which allows access to finite momen- 
tum properties of correlation functions and also allows 
for an accurate calculation of critical exponents. As in 
the derivative expansion, in the NLPA all truncations 
are done at the level of the effective action, a property it 
shares with the derivative expansion. This allows for a 
strict control of the symmetries of the underlying model 
and also allows for extensions of the approach. While 
the present approach includes both terms of arbitrary 
powers in the fields (in the local term), the nonlocal 
terms are restricted up to fourth order in the fields. In 
contrast, in the BMW scheme jllj all vertices have a 
momentum dependence which is however only approxi- 
mately taken into account. The present scheme can easily 
be extended by including for example terms of the type 
p(x) p(y) 2 k(x — y) which would result in momentum de- 
pendent vertices with up to six legs and similar terms of 
higher order in the densities can of course also easily be 
constructed. In D = 3 it might suffice to limit such terms 
only to a small maximal power in p to get converged val- 
ues for the critical exponents. For each such additional 
term a new coupling function must be introduced so the 
nonlocality in the scheme will always be restricted to a 
finite order in the fields. The computational cost of such 
an extension is relatively modest, since one would still 
deal with the flow of a small number of one-parameter 
functions. In contrast, in the BMW the flow must be 
analysed for a two point function which is defined on a 
two-dimensional grid, one dimension each for the depen- 
dence on fields and momenta, which is numerically more 
difficult. 

Even in the the simplest NLPA truncation analysed 
here, the results for both the critical exponents and for 
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TABLE III: Values for the quantity c denned in Eq. (|16[) . from both the symmetric NLPA and the NLPA for the ordered state 
as well as from a perturbative FRG approach (PFRG), the background field scheme (BMW), variational perturbation theory 
(VPT) and Monte Carlo (MC). 



N 


sym. NLPA 


ord. NLPA 


PFRG [9] 


BMW [11] 


VPT [34] 


MC 


1 


1.38 






1.15 


1.07(10) 


1.09(9) [35] 


2 


1.49 


1.60 


1.23 


1.37 


1.27(10) 


1.29(5) [36J 














1.32(2) [37] 


3 


1.59 


1.72 




1.50 


1.43(11) 




4 


1.68 


1.82 




1.63 


1.54(11) 


1.60(10) [35] 


10 


2.02 


2.11 




2.02 







the momentum dependence of the two-point function are 
already surprisingly good and the present scheme offers a 
direct access to both universal and non-universal quanti- 
ties. The present scheme is also certainly useful for more 
complex models where even the local terms are restricted 



to a finite order in the fields 

We thank Pawel Jakubczyk, Andreas Eberlein, and 
Federico Benitez for discussions. This work was sup- 
ported by the DFG research group FOR 723. 
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